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The aim of this mini-course is to give an introduction in Optimal Net- 
works theory. Optimal networks appear as solutions of the following natural 
problem: How to connect a finite set of points in a metric space in an optimal 
way? We cover three most natural types of optimal connection: spanning trees 
(connection without additional road forks), shortest trees and locally shortest 
trees, and minimal fillings. 



1 Introduction: Optimal Connection 

This mini-course was given in the First Yaroslavl Summer School on Discrete and Com- 
putational Geometry in August 2012, organized by International Delaunay Laboratory 
"Discrete and Computational Geometry" of Demidov Yaroslavl State University. We are 
very thankful to the organizers for a possibility to give lectures their and to publish this 
notes, and also for their warm hospitality during the Summer School. The real course 
consisted of three 1 hour lectures, but the division of these notes into sections is inde- 
pendent on the lectures structure. The video of the lectures can be found in the site of 
the Laboratory (http://dcglab.uniyar.ac.ru). The main reference is our books pQ 
and [2] , and the paper [3] for Section [5] 

Our subject is optimal connection problems, a very popular and important kind of 
geometrical optimization problems. We all seek what is better, so optimization problems 
attract specialists during centuries. Geometrical optimization problems related to inves- 
tigation of critical points of geometrical functionals, such as length, volume, energy, etc. 
The main example for us is the length functional, and the corresponding optimization 
problem consists in finding of length minimal connections. 

1.1 Connecting Two Points 

If we have to points A and B in the Euclidean plane R 2 , then, as we know from the 
elementary school, the shortest curve joining A and B is unique and coincides with the 
straight segment AB, so optimal connection problem is trivial in this case. But if we 
change the way of distance measuring and consider, for example, so-called Manhattan 
plane, i.e. the plane R 2 with fixed standard coordinates (a;, y) and the distance function 
Pi{A, B) = \a\ — b\\ + \a 2 ~b 2 \, where A = (aj., a 2 ) and B = (bi, b 2 ), then it is not difficult 
to verify that in this case there are infinitely many shortest curves connecting A and B. 
Namely, if < a\ < b\ and < a 2 < b 2 , then any monotonic curve j(t) = (x(t), y(t)) , 
t £ [0,1], 7(0) = A, 7(1) = B, where functions x(t) and y(t) are monotonic, are the 
shortest, see Figure [T] left. Another new effect that can be observed in this example is 
as follows. In the Euclidean plane a curve such that each its sufficiently small piece is a 
shortest curve joining its ends (so-called locally shortest curve) is a shortest curve itself. 
In the Manhattan plane it is not so. The length of a locally shortest curve having the 
form of the letter II, see Figure [TJ right, can be evidently decreased. 

Similar effects can be observed in the surface of standard sphere S 2 cl 3 . Here the 
shortest curve joining a pair of points is the lesser arc of the great circle (the cross-section 
of the sphere by a plane passing through the origin). Two opposite points are connected 
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Figure 1: The shortest curves connecting a pair of points in Manhattan plane (left), and 
locally shortest but not shortest curve in this plane (right). 

by infinitely many shortest curves, and if points A and B are not opposite, then the 
corresponding great circle is unique and it is partitioned into two arcs, both of them are 
locally shortest, one is the (unique) shortest, but the other one is not. (Really speaking, 
the difference with the Manhattan plane consists in the fact that for the case of the 
sphere any directional derivative of the length of any locally shortest arc with respect to 
its deformation preserving the ends is equal to zero). 

Exercise 1.1 For a pair of points on the surface of the cube describe shortest and locally 
shortest curves. Find out an infinite family of locally shortest curves having pairwise 
distinct lengths. 

1.2 Connecting Many Points: Possible Approaches 

Let us consider general situation, when we are given with a finite set M — {Ai, . . . , A n } 
of points in a metric space {X, p), and we want to connect them in some optimal way in 
the sense of the total length of the connection. We are working under assumption that we 
already know how to connect pairs of points in [X, p), therefore we need just to organize 
the set of shortest curves in appropriate way. There are several natural statements of 
the problem, and we list here the most popular ones. 

1.2.1 No Additional Forks Case: Spanning Trees 

We do not allow additional forks, that is, we can switch between the shortest segments 
at the points from M only. As a result, we obtain a particular case of Graph Theory 
problem about minimal spanning trees in a connected weighted graph. We recall only 
necessary concepts of Graph Theory, the detail can be found, for example in [4] . 

Recall that a [simple) graph can be considered as a pair G = (V, E), consisting of a 
finite set V = {v±, . . . , v n } of vertices and a finite set E = {e\, . . . , e m } of edges, where 
each edge ei is a two-element subset of V. If e = {v, v'}, then we say that v and v' are 
neighboring, edge e joins or connects them, the edge e and each of the vertices v and 
v 1 are incident. The number of vertices neighboring to a vertex v is called the degree 
of v and is denoted by degv. A graph H = (Vh,Eh) is said to be a subgraph of a 
graph G — (Vc, Eg), if Vh C Vg and Ejj C Eg- The subgraph H is called spanning, if 
V H = V g . 

A path 7 in a graph G is a sequence v, x ,e^,Vi 2 . . . , ei k Vi k+1 of its vertices and edges 
such that each edge ei s connects vertices Vi s and v is+1 . We also say that the path 7 
connects the vertices Vi t and Vi k+l which are said to be ending vertices of the path. A 
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Figure 2: Minimal spanning tree (left), shortest tree (center), and minimal filling, con- 
necting the vertex set of regular triangle in Euclidean plane. 

path is said to be cyclic, if its ending vertices coincide with each other. A cyclic path 
with pairwise distinct edges is referred as a simple cycle. A graph without simple cycles 
is said to be acyclic. A graph is said to be connected, if any its two vertices can be 
connected by a path. An acyclic connected graph is called a tree. 

If we are given with a function lj: E — > 1 on the edge set of a graph G, then 
the pair (G,lu) is referred as a weighted graph. For any subgraph H — (Vh,Eh) of a 
weighted graph (G = (Vg, Eh),uj) the value uj(H) = J2eeE H CJ ( e ) ^ s caue d the weight of 

H. Similarly, for any path 7 = 1^, , Vi 2 . . . , ei k Vi k+1 the value 1^(7) ~ X) s =i a; ( e * s ) ^ s 
called the weight ofj. 

For a weighted connected graph (G = {Vg, Eg),<jj) with positive weight function u>, 
a spanning connected subgraph of minimal possible weight is called minimal spanning 
tree. The positivity of u) implies that such subgraph is acyclic, i.e. it is a tree indeed. 
The weight of any minimal spanning tree for (G,ui) is denoted by mst(G, oj). 

Optimal connection problem without additional forks can be considered as minimal 
spanning tree problem for a special graph. Let M — {A\, . . . , A n } be a finite set of points 
in a metric space (X, p) as above. Consider the complete graph K(M) with vertex set M 
and edge set consisting of all two-element subsets of M . In other words, any two vertices 
Ai and Aj are connected by an edge in K(M). By AiAj we denote the corresponding 
edge. The number of edges in K(M) is, evidently, n(n — l)/2. We define the positive 
weight function uj p (AiAj) — p(Ai,Aj). Then any minimal spanning tree T in K(M) 
can be considered as a set of shortest curves in (X, p) joining corresponding points and 
forming a network in X connecting M without additional forks in an optimal way, i.e. 
with the least possible length. Such a network is called a minimal spanning tree for M in 
(X,p). Its total weight w p (T) is called length and is denoted by mstx(-^0- -"- n Section [2] 
we speak about minimal spanning trees in more details. 

I. 2.2 Shortest tree: Fermat Steiner Problem 

But already P. Fermat and C. F. Gauss understood that additional forks can be profitable, 
i.e. can give an opportunity to decrease the length of optimal connection. For example, see 
Figure^ if we consider the vertex set M = {Ai, Ai, A3} of a regular triangle with side 1 
in the Euclidean plane, then the corresponding graph K(M) consists of three edges of the 
same weight 1 and each minimal spanning tree consists of two edges, so mst R 2(M) = 2. 
But if we add the center T of the triangle and consider the network consisting of three 
straight segments A{I ', A 2 T, A 3 T, then its length is equal to 3|-^ = \/3 < 2, so it is 
shorter than the minimal spanning tree. 

This reasoning leads to the following general definition. Let M — {Ai, . . . , A n } be 
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a finite set of points in a metric space (X, p) as above. Consider a larger finite set N, 
M C N C X, and a minimal spanning tree for N in X. Then this tree contains M 
as a subset of its vertex set N, but also may contain some other additional vertices- 
forks. Such additional vertices are referred as Steiner points. Further, we define a value 
$mtx(M) = infpf-.McNcx mstx(N) and call it the length of shortest tree connecting M or 
of Steiner minimal tree for M . If this infimum attains at some set N, then each minimal 
spanning tree for this N is called a shortest tree or a Steiner minimal tree connecting M. 
Famous Steiner problem is the problem of finding a shortest tree for a given finite subset 
of a metric space. We will speak about Steiner problem in more details in Section[3J The 
shortest tree for the vertex set of a regular triangle in the Euclidean plane is depicted in 
Figure [H 



1.2.3 Minimizing over Different Ambient Spaces: Minimal Fillings 

Shortest trees give the least possible length of connecting network for a given finite set in 
a fixed ambient space. But sometimes it s possible to decrease the length of connection 
by choosing another ambient space. Let M — {Ai, . . . , A n } be a finite set of points 
in a metric space {X, p) as above, and consider M as a finite metric space with the 
distance function pm obtained as the restriction of the distance function p. Consider 
an isometric embedding ip: (M,pm) — > (Y, py) of this finite metric space (M,pm) into 
a (compact) metric space (Y, py) and consider the value smty {tp(MJ\. It could be less 
than smtx(M). For example, the vertex set of the regular triangle with side 1 can be 
embedded into Manhattan plane as the set {(-1/2, 0), (0, 1/2), (1/2, 0)}, see Figure [U 
Than the unique additional vertex of the shortest tree is the origin and the length of 
the tree is 3/2 < \/3. So, for a finite metric space Ai = (M,pm), consider the value 
mf(A^) = inf^ smty- ((p(M)) which is referred as weight of minimal filling of the finite 
metric space Ai. Minimal fillings can be naturally defined in terms of weighted graphs 
and can be considered as a generalization of Gromov's concept of minimal fillings for 
Riemannian manifolds. We speak about them in more details in Section [5l 



2 Minimal Spanning Trees 

In this section we discuss minimal spanning trees construction in more details. As we have 
already mentioned above, in this case the problem can be stated in terms of Graph Theory 
for an arbitrary connected weighted graph. But geometrical interpretation permits to 
speed up the algorithms of Graph Theory. 



2.1 General Case: Graph Theory Approach 

We start with the Graph Theory problem of finding a minimal spanning tree in a con- 
nected weighted graph. It is not difficult to verify that direct enumeration of all possible 
spanning subtrees of a connected graph leads to an exponential algorithm. 

To see that, recall well-known Kirchhoff theorem counting the number of span- 
ning subtrees. If G = (V, E) is a connected graph with enumerated vertex set V — 
{vi, . . . , «„}, then its Kirchhoff matrix is defined as (n x n)-matrix Bq — (6y) with 
elements 

deg-^i if * = i, 
-1 if {vi, Vj} G E, 

otherwise. 



Then the following result based on elementary Graph Theory and Binet-Cauchy formula 
for determinant calculation is valid, see proof, for example in [3]. 
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Theorem 2.1 (Kirchhoff) For a connected graph G with n > 2 vertices, the number 
of spanning subtrees is equal to the algebraic complement of any element of the Kirchhoff 
matrix Bq . 

Example. Let G = K n be the complete graph with n vertices. Than its Kirchhoff 
matrix has the following form: 



B 



K„ 



I n- 1 -1 -1 

-1 71-1-1 
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-1 -1 -1 •■■ 

The algebraic complement of the element b nn is equal to 
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where the first equality is obtained by change of the first row by the sum of all the rows, 
and the second equality is obtained by change of the ith row, i > 2, by the sum of it 
with the first row. 



Corollary 2.2 The complete graph with n vertices contains n n 2 spanning trees. 

Remark. Notice that this result is equivalent to Cayley Theorem saying that the total 
number of trees with n enumerated vertices is equal to n n . 

But it is a surprising fact, that there exist polynomial algorithms constructing minimal 
spanning trees. Several such algorithms were discovered in 1960s. We tell about Kruskal's 
algorithm. Similar Prim's algorithm can be found in [4]. 

So, we are given with a connected weighted graph (G — (V,E),ui) with positive 
weight function oj. At the initial step of Kruckal algorithm we construct the graph 
To = (V, 0) and put E = E. If the graph T,_i and the non-empty set C E, 

i < n — 1, have been already constructed, then we choose in Ei-i an edge ej of least 
possible weight and construct a new graph T = T%-\ U ej and also a new set Ei = {e G 
E | e $Ti and T U e is acyclic}. Algorithm stops when the graph T n _i is constructed^. 

Theorem 2.3 (Kruskal) Under the above notations, the graph T„_i can be constructed 
for any connected weighted graph Q = (G,uj), and moreover, T„_i is a minimal spanning 
tree in Q . 

Proof. The set Ei is non-empty for all < i < n — 2, because the corresponding 
subgraphs Ti are not connected (the graph Ti has n vertices and i edges), therefore all 
the graphs T\, . . . , T„_i can be constructed. Further, all these graphs are acyclic due to 
the construction, and T„_i has n vertices and n — 1 edges, so it is a tree. 

To finish the proof it remains to show that the spanning tree T„_i C G is minimal. 
Since the graph Q has a finite number of spanning trees, a minimal spanning tree does 
exist. Let T be a minimal spanning tree. We show that it can be reconstructed to the 
tree T„__i without changing the total weight, so T„_i is also a minimal spanning tree. 

1 Here the operation of adding an edge e to a graph G = (V, E) can be formally defined as follows: 
G U e = (V, E U {e}). Similarly, G\e=(V,E\ {e}). 
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To do this, recall that the edges of the tree T„_i are enumerated in accordance with 
the work of the algorithm. Denote them by e%, . . . , e n -i as above, and assume that efc 
is the first one that does not belong to T . The graph T U efc contains a unique cycle 
c D et- This cycle c also contains an edge e not belonging to T n _i (otherwise c C T n _i, 
a contradiction). Consider the graph T" = T U efc \ e. It is evidently a spanning tree in 
G, and therefore its weight is not less than the weight of the minimal spanning tree T, 
hence 

w(T') = w(T) + w(e fc ) - w(e) > w(T), 

and thus, w(efc) > w(e). 

On the other hand, all the edges e%, . . . , efc_i belongs to T by our assumption. There- 
fore, the graph Xfc_iUe is a subgraph of T and is acyclic, in particular. Hence, e € -Efc-i so 
as efc e Ek-i- But the algorithm has chosen efc, hence w(efc) < ui(e). Thus, u)(ek) — u(e), 
and so w(T') = lj(T), and therefore T' is a minimal spanning tree in C?. But now T" 
contains the edges ei, . . . , efc from T n _i. Repeating this procedure we reconstruct T to 
T„_i in the class of minimal spanning trees. Theorem is proved. 

Remark. For a connected weighted graph with n vertices and m edges the complexity 
of the Kruskal's algorithm can be naturally estimated as mn ~ n 3 . The estimation 
can be improved to mlogm ~ n 2 \ogn. The fastest non-randomized comparison-based 
algorithm with known complexity belongs to Bernard Chazelle [5]- It turns out that if 
the weight function is geometrical, then the algorithms can be improved. 

2.2 Euclidean Case: Geometrical Approach 

Now assume that M is a finite subset of the Euclidean plane K 2 . It turns out that a 
minimal spanning tree for M in M 2 can be constructed faster than the one for an abstract 
complete graph with n = \M\ vertices by means of some geometrical reasonings. To do 
that we need to construct so called Voronoi partition of the plane, corresponding to M, 
and the Delaunay graph on M. It turns out that any minimal spanning tree for M in 
R 2 is a subgraph of the Delaunay graph, see Figure El and the number of edges in this 
graph is linear with respect to n, so the standard Kruskal's algorithm applied to it gives 
the complexity nlogn instead of n 2 logn for the complete graph with n vertices. 

Let us pass to details. Let M = {A±, . . . ,A n } c M 2 be a finite subset of the plain. 
The Voronoi cell of the point Ai is defined as 

Vor M (^) = {x e M 2 | ||x - A{ || < \\x - AA\ for all j} . 

The Voronoi cell for Aj is a convex polygonal domain which is equal to the intersection 
of the closed half-planes restricted by the perpendicular bisectors of the segments Aj,Aj, 
j =/= i. It is easy to verify, that the intersection of any two Voronoi cells has no interior 
points and that Uj VorM(A^) = M 2 . This partition of the plane is referred as Voronoi 
partition or Voronoi diagram. Two cells VovM(Ai) and Vorjvf (Aj) are said to be adjacent, 
if there intersection contains a straight segment. The Delaunay graph D(M) is defined 
as the dual planar graph to the Voronoi diagram. More precisely, the vertex set of D(M) 
is M, and to vertices Ai and Aj are connected by an edge, if and only if their Voronoi 
cells Voi M(Ai) and Vor M(Aj) are adjacent. The edges of the Delaunay graph are the 
corresponding straight segments. 

It is easy to verify, that if the set M is generic in the sense that no three points lie 
at a common straight line and no four points lie at a common circle, then the Delaunay 
graph D(M) is a triangulation, i.e. its bounded faces are triangles. In general case some 
bounded faces could be inscribed polygons. Anyway, the number of edges of the graph 
D{M) does not exceed 3n. It remains to prove the following key Lemma. 
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Figure 3: Voronoi diagram (left) and Delaunay graph together with a minimal spanning 
tree (right) for a point set in the plane. 

Lemma 2.4 Any minimal spanning tree for M C R 2 is a subgraph of the Delaunay 
graph D(M). 

Proof. Let e = AiAj be an edge of a minimal spanning tree T for M. We have to show 
that the Voronoi cells VoiM(Ai) and Yotm(Aj) are adjacent. The graph T\e consists 
of two connected components, and this partition generates a partition of the set M into 
two subsets, say Mi and M%. Assume that Ai <G M\ and Aj 6 M%. The minimality of 
the spanning tree T implies that ||A, A,-| is equal to the distance between the sets Mi 
and M2, where stands for the distance between Ai and Ak- 

By u we denote the middle point of the straight segment AiAj, and let Ak be another 
point fromM. Assume that Ak G M 2 . Due to the previous remark, > H^A,-]], 

therefore 

\\u,A k \\ > \\Ai,A k \\ - \\Ai,u\\ > \\AiAj\\ - \\Ai,u\\ = H^, A,-||/2 = \\u,Ai\\ = \\u,Aj\\. 

On the other hand, if \\u, Ak\\ = \\u, Ai\\, then we have equalities in both above inequali- 
ties. The first one means that u lies at the straight segment AiAk, hence Ak lies at the ray 
AiU. The second equality implies HAiAfcH = ||Ajj4j||, and so Ak = Aj, a contradiction. 
Thus, \\u, Ak\\ > \\u,Ai\\, that is u does not belong to the cell Vorjvf(Afe) for k £ {i, j}. 
Thus, u € Vor M (Ai) n Vor A/ (A :) ), 

Since the inequality proved is strict, the same arguments remain valid for points lying 
close to u on the perpendicular bisector to the segment A^Aj. Therefore, the intersection 
of the Voronoi cells VorM^i) and VorM (Aj) contains a straight segment, that is the 
cells are adjacent. Lemma is proved. 

Remark. The previous arguments work in any dimension. But the trouble is that 
starting from the dimension 3 the number of edges in the Delaunay graph need not be 
linear on the number of its vertices. 

Exercise 2.5 Verify that the same arguments can be applied to minimal spanning trees 
for a finite subset of M™ . 

Exercise 2.6 Give an example of a finite subset Mcl 3 such that the Delaunay graph 
D(M) coincides with the complete graph K(M). 

Problem 2.7 In what metric spaces similar geometrical approach also works? It defi- 
nitely works for planar polygons with intrinsic metric, see Jfjjj. 
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Figure 4: Torricelli-Simpson construction, the case Ai < 120° (left), and the case A3 > 
120° (right). 

3 Steiner Trees and Locally Minimal Networks 

In this section we speak about shortest trees and locally shortest networks in more de- 
tails. Besides necessary definitions we discuss local structure theorems, Melzak algorithm 
constructing locally minimal trees in the plane, global results concerning locally minimal 
binary trees in the plane (so called twisting number theory) and the particular case, lo- 
cally minimal binary trees with convex boundaries (language of triangular tilings) . The 
details concerning twisting number and tiling realization theory can be found in [2] or [1] , 
and also in [7]. 

3.1 Fermat Problem 

The idea that additional forks can help to decrease the length of a connecting network 
had been already clear to P. Fermat and his students. It seems that Fermat was the 
first, who stated the following optimization problem: for given three points A\, A2, and 
^3 in the plane hnd a point X minimizing the sum of distances from the points Ai, 
i.e. minimize the function F(X) — For the case when all the angles of the 

triangle A1A2A3 are less than or equal to 120° the solution was found by E. Torricelli 
and later by R. Simpson. The construction of Torricelli is as follows, see Figure 2] 

On the sides of the triangle A1A2A3 construct equilateral triangles AiAjA' k , {i,j, k} = 
{1, 2, 3}, such that they intersect the initial triangle only by the common sides. Then, as 
Torricelli proved, the circumscribing circles of these three triangles intersect in a point 
referred as Torricelli point T of the triangle A1A2A3. If all the angles Ai are less than 
or equal to 120°, then T lies in the triangle A1A2A3 and gives the unique solution to 
the Fermat problem^ Later Simpson proved that the straight segments AiA\ also pass 
through the Torricelli point, and the lengths of all these three segments are equal to 
F(T). If one of the angles, say A3, is more than 120°, then the Torricelli point is located 
outside the triangle and can not be the solution to Fermat problem. In this case the 
solution is X = A3. 

2 An elementary proof can be obtained by rotation R of a copy of the triangle around its vertex, say 
Ai, by 60° and considering the polygonal line L joining A2, X, image R(X) of X under the rotation, 
and R,(As). The length of L is equal to F(X), and minimal value of F(X) corresponds to the location 
of the X such that L is a straight segment. 
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Remark. So we see, that shortest tree for a triangle in the plane consists of straight 
segments meeting at the vertices by angles more than or equal to 120°. It turns out, that 
this 120° -property remains valid in much more general situation. 

3.2 Local Structure Theorem and Locally Minimal Networks 

Let M = {A 1} . . . , A n } be a finite subset of Euclidean space M. N , and T is a Steiner tree 
connecting M, Recall that we defined shortest trees as abstract graphs with vertex set in 
the ambient metric space. In the case of M. N it is natural to model edges of such graph as 
straight segments joining corresponding points in the space. The configuration obtained 
is referred as a geometrical realization of the corresponding graph. Below, speaking about 
shortest trees in 1^ we will usually mean their geometrical realizations. The local 
structure of a shortest tree (more exactly of a geometrical realization of the tree) can be 
easily described. 

Theorem 3.1 (Local Structure) Let V be a shortest tree connecting a finite subset 
M = {Ai, . . . , A n } mR N . Then 

1. all edges ofT are straight segments] 

2. any vertex v £ T of degree 1 belongs to M; 

3. any two neighboring edges ofT meet in common vertex by angle more than or equal 
to 120°; 

4- if the degree of a vertex v is equal to 2 and v (jL M, then the edges meet at v by 
180° angle. 

Corollary 3.2 Let T be a shortest tree connecting a finite subset M — {A\, . . . ,A n } in 
M. N . Then the degree of any its vertex is at most 3, and if the degree of a vertex v equals 
to 3, then the edges meet at v by angles equal to 120°. 

Example. Let M be the vertex set of regular tetrahedron A in R 3 . Then the network 
consisting of four straight segments joining the vertices of the tetrahedra with its center 
O is not a shortest network. Indeed, since degO = 4, then the angles between the 
edges meeting at O are less than 120°. The set M is connected by three different (but 
isometrical) shortest networks, each of which has two additional vertices of degree 3, see 
Figure [5] 

Theorem 13 . 1 1 can be just "word- by- word" extended to the case of Riemannian mani- 
folds (we only need to change straight segments by geodesic segments) [T] and even to the 
case of Alexandrov spaces with bounded curvature. The case of normed spaces turned 
out to be more complicated (some general results can be found in [2]). 

A connected graph V in M. N (in a Remannian manifold) whose vertex set contains 
a finite subset M C is called a locally minimal network connecting M or with the 
boundary dT = M, if it satisfies Conditions (l)-(4) from Theorem 13.11 In the case of 
complete Remannian manifolds such graphs are minimal "in small," i.e. the following 
result holds, see [T]. 

Theorem 3.3 (Minimality "in small") Let T be a locally minimal network connect- 
ing a finite subset M of a complete Riemannian manifold W . Then each point P G T 
possesses a neighborhood U in W , such that the network TnU is a shortest network with 
the boundary (dT f] U) U (r n dU). 

Remark. In the case of normed spaces Theorem 13. 31 is not valid even for two-point sets, 
see example in Figure [TJ 
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Figure 5: Non-shortest tree (left) and one of the shortest trees (right) for the vertex set 
of a regular tetrahedron. 
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Figure 6: Tree G is partitioned into 4 non-degenerate components by cutting at vertices 
of degree 2. 

3.3 Melzak Algorithm and Steiner Problem Complexity 

Let us return back to the case of Euclidean plane. It turns out that in this case the 
Torricelli-Simpson construction can be generalized to a geometrical algorithm, that either 
constructs a locally minimal tree of a given structure for a given boundary set, or reports 
that such a tree does not exist. This algorithm was discovered by Z. Melzak [9] and 
improved by F. Hwang [5]. 

Assume that we are given with a tree G whose vertex degrees are at most 3, a finite 
subset M of the plane, and a bijection ip : dG — > M, where dG is the set of all vertices 
from G of degrees 1 and 2. To start with, partition the tree G into the union of so-called 
non-degenerate components Gi by cutting the tree at each its vertex of degree 2, see 
Figure O To construct locally minimal network T of type G spanning M in accordance 
with <p it suffices to construct each its component Yi of type Gi on the corresponding 
boundary Mj = ip(dGi), where dGi = dG n G;, in accordance with (fi — ip\dGi and to 
verify the angles between the edges of the components at the vertices of degree 2. All 
these angles must be more than or equal to 120°, see Figured] 

Now we pass to the case of one non-degenerate component, i.e. we assume that G 
has no vertices of degree 2 and that dG consists of all the vertices of degree 1. Such 
trees are referred as binary. If \dG\ = 2, then the corresponding locally minimal tree Y 
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Figure 7: A step of forward trace of Melzak algorithm. 



is a straight segment. Otherwise, it is easy to verify that each such tree G contains so- 
called moustaches, i.e. a pair of vertices of degree 1 neighboring with a common vertex 
of degree 3. Fix such moustaches {x, x'} C dG, by y denote their common vertex of 
degree 3, and make a forward step of Melzak algorithm, see Figure [7J that reduces the 
number of boundary vertices by 1. Namely, we reconstruct the tree G by deleting the 
vertices x and x' together with the edges xy and x'y and adding y to the boundary of new 
binary tree; reconstruct the set M by deleting the points ip(x) and <p(x') and adding a 
new point A xx i which is the third vertex of a regular triangle constructed on the straight 
segment ip(x)(p(x') in the plane; and reconstruct the mapping tp putting ip(y) = A xx >. 
Notice that the point A xx i can be constructed in two ways, because there are two such 
regular triangles. Thus, if the number of boundary vertices in the resulting tree is more 
than 2, then we can repeat the procedure described above. And if it becomes 2, then we 
can construct the corresponding locally minimal tree — the straight segment. Here the 
forward trace of Melzak algorithm stops. Now we have to reconstruct the initial tree, if 
possible. 

Thus, we have a straight segment /Cl 2 realizing locally minimal tree with unique 
edge ab, and at least one of its ending points has the form A xx >, where x and x' are the 
boundary vertices of the binary tree G from the previous step, neighboring with their 
common vertex of degree 3. Let this common vertex be a, that is a corresponds to A xx > . 
We reconstruct G by adding edges ax and ax' . Then we restore the points f(x) and 
ip(x') in the plane together with the regular triangle (p(x)ip(x')A xx > , circumscribe the 
circle S 1 around it and consider the intersection of S 1 with the segment /, see Figure [8] 
If it does not contains a point lying at the smaller arc of S 1 restricted by ip(x) and ip(x'), 
then the tree G can not be reconstructed and we have to pass to another realization of 
the forward trace of the algorithm. Otherwise we put ip{a) be equal to this point. The 
straight segments (p(x)(p(a) and <p(x')(p(a) meet at (p(a) by 120° and together with the 
subsegment tp(a)(p(b) form a locally minimal binary tree T of type G with tree boundary 
vertices. We repeat this procedure until we either reconstruct the tree of type G, or 
verify all possible realizations of the forward trace and conclude that the tree of type G 
does not exists. 

The Melzak algorithm described above contains an exponential number of possibili- 
ties of its forward trace realization, due to two possible locations of each regular triangle 
constructed by the algorithm. This complexity can be reduced by means of modification 
suggested by F. Hwang [8]. He showed that considering a bit more complicated config- 
urations of boundary points (four points corresponding to "neighboring moustaches" or 
three points corresponding to moustaches and "neighboring" degree-1 vertex) one always 
can understand which regular triangle must be chosen, see details in [8]. 

But unfortunately even a linear time realization of Melzak algorithm does not lead to 
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Figure 8: A step of back trace of Melzak algorithm. 




Figure 9: Three binary trees with 4 vertices of degree 1. 

a polynomial algorithm of a shortest tree finding. The reason is a huge number of possible 
structures of the tree G with \dG\ — n together with also exponential number of different 
mappings ip: dG — > M for fixed dG and M . Even for binary trees we have 3 possibilities 
for n — 4, see Figure [9j and 15 possibilities for n = 5 (notice that the corresponding 
binary trees are isomorphic as graphs). For n = 6 we have two non- isomorphic binary 
trees and the number of possibilities becomes 90. It can be shown that the total number 
of possibilities can be estimated by Catalan number and grough exponentially. 

So, to obtain an efficient algorithms, we have to find some a priori restrictions on pos- 
sible structures of minimal networks. In the next subsection we tell about the restrictions 
generated by geometry of boundary sets. 

3.4 Boundaries Geometry and Networks Topology- 
Here we review our results from jTD] and [7]. The goal is to find some restriction on 
the structure of locally minimal binary trees spanning a given boundary in the plane in 
terms of geometry of the boundary set. To do this we need to choose or to introduce 
some characteristics of the network structure and of the boundary geometry. 

As a characteristic of the geometry of a boundary set M we take the number of 
convexity levels c(M). Recall the definition. Let M be a finite non-empty subset of 
the plane. Take the convex hull chM of M and assign the points from M lying at the 
boundary of the polygon chM to the first convexity level M^- 1 ' of M. If the set M\M^ 
is not empty, then define the second convexity level of M to be equal to the first 

convexity level of M \ M^ x \ and so on. As a result, we obtain the partition of the set 
M into its convexity levels, and by c(M) we denote the total number of this levels, see 
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Figure 10: A boundary set with 4 convexity levels. 

Figure HD1 

Now let us pass to definition of a characteristic describing the "complexity" of planar 
binary trees. Assume that we are given with a planar binary tree T, and let the orientation 
of the plane be fixed. For any its two edges, say e\ and e2, we consider the unique path 
7 in r starting at e\ and finishing at e%- All interior vertices of 7 are the vertices of Y 
having degree 3. Let us walk from e\ to ei along 7. Then at each interior vertex of 7 
we make either left, or right turn in Y. Define the value tw(ei,e2) to be equal to the 
difference between the numbers of left and right turns we have made. In other words, 
assign to an interior vertex of 7 the label r = ±1, where +1 corresponds to left turns and 
— 1 to right turns. Then tw(ei, e-i) is the sum of these values, see Figure [TT] Notice that 
tw(ei,e2) = — tw(e2,ei). At last, we put twT = maxtw(ej, ej), where the maximum is 
taken over all ordered pairs of edges of Y. 

If the tree Y is locally minimal, then the twisting number between any pair of its 
edges has a simple geometrical interpretation, see Figure 1111 Namely, since the angles 
between any neighboring edges are equal to 27r/3, then tw(ej,ej) is equal to the total 
angle which the oriented edge rotates by passing from to ej, divided by it/3. 

It turns out, that the twisting number of a locally minimal binary tree with a given 
boundary is restricted from above by a linear function on the number of convexity levels 
of the boundary. Namely, the following result holds. 

Theorem 3.4 Let Y be a locally minimal binary tree connecting the boundary set M that 
coincides with the set of vertices of degree 1 from Y. Then 

twY < 12 (c(M) - l) + 5. 

The important particular case c{M) = 1 corresponds to the vertex sets of convex 
polygons. Such boundaries are referred as convex. 

Theorem 3.5 Let Y be a locally minimal binary tree with a convex boundary. Then 
twT < 5. Conversely, any planar binary tree Y with twT < 5 is planar equivalent to a 
locally minimal binary tree with a convex boundary. 

Notice that the direct statement of Theorem 13.51 is rather easy to prove (it follows 
from the geometrical interpretation of the twisting number, easy remark that twT al- 
ways attains at boundary edges, and the monotony of convex polygonal lines). But the 
converse statement is quite non-trivial. The proof obtained in [7] is based on the com- 
plete description of binary trees with twisting number at most five, obtained in terms of 
so-called triangular tilings that will be discussed in the next subsection. 
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Figure 11: Left and right turns in planar binary tree and its twisting number (left), and 
twisting number of locally minimal binary tree. 



Figure 12: A diagonal triangulation and corresponding planar binary tree. 



Problem 3.6 Estimate the number of binary trees structures with n vertices of degree 1 
and twisting number at most k. It is more or less clear that the number is exponential 
on n even for k — 5, but it is interesting to obtain an exact asymptotic. 

3.5 Triangular Tilings and their Applications 

It turns out that the description of planar binary trees with twisting number at most five 
can be effectively done in the language of planar triangulations of a special type which 
are referred as triangular tilings. 

The correspondence between diagonal triangulations of planar convex polygons and 
planar binary trees is well-known: the planar dual graph of such triangulation is a binary 
tree, see Figure [12J and each binary tree can be obtained in such a way. Here the vertices 
of the dual graph are centers of the triangles of the triangulation (medians intersection 
point) and middle points of the sides of the polygon; and edges are straight segments 
joining either the middle of a side with the center of the same triangle, or two centers of 
the triangles having a common side. 

In the context of locally minimal binary trees, the most effective way to represent 
the diagonal triangulations is to draw them consisting of regular triangles. Such special 
triangulations are referred as triangular tilings. The main advantage of the tilings is that 
the dual binary tree constructed as described above is a locally minimal binary tree with 
the corresponding boundary. Therefore, tilings "feel the geometry" of locally minimal 
binary trees and turns out to be very useful in the description of such trees with small 
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Figure 13: A triangular tiling with its dual binary tree (left), its outer cells (middle) and 
inner cells (right). 




Figure 14: Un-paired and paired growths of a tiling (left), growths that should be 
deleted to get a skeleton (middle) , corresponding decomposition into skeleton and growths 
(right). 

twisting numbers. 

The main difficulty in constructing a triangulation consisting of regular triangles for 
a given binary tree is that the resulting polygon can overlap itself. An example of such 
overlapping can be easily constructed from a binary tree T corresponding to the diagonal 
triangulation of a convex n-gon, n > 6, all whose diagonals are incident to a common 
vertex. But the twisting number of such T is also at least 6. The following result is 
proved in [7]. 

Theorem 3.7 The triangular tiling corresponding to any planar binary tree with twisting 
number less than or equal to five has no self-intersections. 

Theorem 13.71 gives an opportunity to reduce the description of the planar binary trees 
with twisting number at most five to the description of the corresponding triangular 
tilings. 

To describe all the triangular tilings whose dual binary trees have the twisting number 
at most five, we decompose each such tiling into elementary "breaks" . The triangles of 
the tiling are referred as cells. A cell of a tiling T is said to be outer, if two its sides lie 
at the boundary of T considered as planar polygon. Further, a cell is said to be inner, 
if no one of its sides lies at the boundary, see Figure [13] An outer cell adjacent to (i.e. 
intersecting with by a common side) an inner cell is referred as a growth of T. A tiling 
can contains as un-paired growths, so as paired growths, see Figure HU 

For each inner cell we delete exactly one growth adjacent to it, providing such growths 
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Figure 15: Branching points (left), linear parts (middle), and code of a skeleton (right). 
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Figure 16: All possible codes of skeletons whose dual binary trees have twisting number 
at most 5. 

exist. As a result, we obtain a decomposition of the initial tiling into its growths and its 
skeleton (a tiling without growths). Notice, that such a decomposition is not unique. 

It turns out that the skeletons of the tilings whose dual binary trees have twisting 
number at most five can be described easily. Also, the possible location of growthes in 
such tilings on their skeletons also can be described. The details can be found in [7] or [1]. 
Here we only formulate the skeletons describing Theorem and include several examples 
of its application. 

Inner cells of a skeleton S are organized into so-called branching points, see Figure [TBI 
After the branching points deleting, the skeleton is partitioned into linear parts. Each 
linear part contains at most one outer cell. Construct a graph C(S) referred as the code 
of the skeleton S as follows: the vertex set of C(S) is the set of its branching points 
and of the outer cells of its linear parts. The edges correspond to the linear parts, see 
Figure [151 

The following result is proved in [7]. 

Theorem 3.8 Consider all skeletons whose dual graphs twisting numbers are at most 5 
and for each of these skeletons construct its code. Then, up to planar equivalence, we 
obtain all planar graphs with at most 6 vertices of degree 1 and without vertices of degree 
2. In particular, every such skeleton contains at most 4 branching points and at most 9 
linear parts. 

All possible codes of such skeletons are depicted in Figure [Tf)J 

This description of skeletons and corresponding tilings obtained in [7], was applied 
to the proof of inverse (non-trivial) statement of Theorem 13.51 In some sense, the proof 
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Figure 17: Two infinite families of locally minimal binary trees connecting vertex sets of 
regular polygons, existing for any regular n-gon (left) and for 3k + 6-gons (middle), and 
the unique finite family, existing for 24-, 30-, 36- and 42-gons only (right). 

obtained in [7] is constructive: for each tiling under consideration a corresponding locally 
minimal binary tree with a convex boundary is constructed. 

Another application is a description of all possible binary trees of the skeleton type 
that can be realized as locally minimal binary trees connecting the vertex set of a regular 
polygon. It turns out, see details in pQ, that there are 2 infinite families of such trees and 
1 finite family. The representatives of these networks together with the corresponding 
skeletons are shown in Figure [T7J 

4 Steiner Ratio 

As we have already discussed in the previous Section, the problem of finding a shortest 
tree connecting a given boundary set is exponential even in two-dimensional Euclidean 
plane. On the other hand, in practice it is necessary to solve transportation problems of 
this kind for several thousands boundary points many times a day. Therefore, in practice 
some heuristical algorithms are used. One of the most popular heuristics for a shortest 
tree is corresponding minimal spanning tree. But using such approximate solutions 
instead of exact one it is important to know the value of possible error appearing under 
the approximation. The Steiner ratio of a metric space is just the measure of maximal 
possible relative error for the approximation of a shortest tree by the corresponding 
minimal spanning tree. 

4.1 Steiner Ratio of a Metric Space 

Let M be a finite subset of a metric space (X,p), and assume that \M\ > 2. We put 
sr M = smt(M) / mst(M). Evidently, sr M < 1. The next statement is also easy to prove. 

Assertion 4.1 For any metric space (X,p) and any its finite subset M C X, \M\ > 2, 
the inequality srM > 1/2 is valid. 

Proof. Let G be a Steiner tree connecting M, Consider an arbitrary embedding of the 
graph G into the plane, walk around G in the plane and list consecutive paths forming 
this tour and joining consecutive boundary vertices from M. The length of each such 
path 7pq joining boundary vertices PQ, i.e. the sum of the lengthes of its edges, is more 
than or equal to the distance p(P, Q), due to the triangle inequality. Consider the cyclic 
path in the complete graph with vertex set M consisting of edges formed by the pairs 
of consecutive vertices from the tour, and let T be a spanning tree on M contained in 
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this path. It is clear, that p(T) < J2(p q) p(Jpq)> where the summation is taken over 
all the pairs of consecutive vertices of the tour. On the other hand, each edge of the 
tree G belongs to exactly two such paths, hence J2(pq) p(1pq) = 2p(G). So, we have 
sr(M) > p(G)/p(T) > 1/2. The Assertion is proved. 

The value sr(M) is the relative error appearing under approximation of the length of 
a shortest tree for a given set M by the length of a minimal spanning tree. The Steiner 
ratio of a metric space (X,p) is defined as the value sr(X) — infMcxsr(M), where the 
infimum is taken over all finite subsets M, \M\ > 2 of the metric space X. So, the Steiner 
ratio of X is the value of the relative error in the worse possible case. 

Corollary 4.2 For arbitrary metric space (X,p) the inequality 1/2 < sr(X) < 1 is valid. 

Exercise 4.3 Verify, that for any r G [1/2,1] there exists a metric space (X,p) with 
sr(X ) = r , see corresponding examples in f^. 

Sometimes, it is convenient to consider so-called Steiner ratios sr n (X) of degree n, 
where n > 2 is an integer, which are defined as follows: sr n (A) = mf M cx,|M|<n sr(M). 
Evidently, Si2{X) = 1. It is also clear that sr(X) = inf n sr n (X). 

Steiner ratio was firstly defined for the Euclidean plane in [TTj , and during the follow- 
ing years the problem of Steiner ratio calculation is one of the most attractive, interesting 
and difficult problems in geometrical optimization. A short review can be found in 2 
and in [12] . One of the most famous stories here is connected with several attempts to 
prove so-called Gilbert-Pollack Conjecture, see [11], saying that sr(R 2 , p-i) = V3/2, where 
pi stands for the Euclidean metric, and hence sr(R 2 , p-i) is attained at the vertex set of a 
regular triangle, see Figure [2] In 1990s D. Z. Du and F. K. Hwang announced that they 
proved the Steiner Ratio Gilbert-Pollak Conjecture [13] , and their proof was published in 
Algorithmica [14]. In spite of the appealing ideas of the paper, the questions concerning 
the proof appeared just after the publication, because the text did not appear formal. 
And about 2003-2005 it becomes clear that the gaps in the D. Z. Du and F. K. Hwang 
work are too deep and can not be repaired, see detail in [15] . 

4.2 Steiner Ratio of Small Degrees for Euclidean Plane 

Gilbert and Pollack calculated sr3(R 2 ,p2) in their paper [11]. We include their proof 
here. 

Since the Steiner ratio of a regular triangle is equal to V3/2, then sr 3 (R 2 , p 2 ) < V3/2, 
so we just need to prove the opposite inequality. To do this, consider a triangle ABC 
in the plane. If one of its angles is more than or equal to 120°, then the shortest tree 
coincides with minimal spanning tree, so in this case sr( ABC) = 1. So it suffices to 
consider the case when all the angles of the triangle are less than 120°. 

Let S be the Torricelli point of the triangle ABC. Show firstly that \AS\ < \BS\, if 
and only if \BC\ > \AC\, i.e. the shortest edge of the Steiner minimal tree lies opposite 
with the longest side of the triangle. The proof is shown in Figure [TBI left- Indeed, if 
\BS\ < \ AS\, then we take the point B' G [S,B] with \SB'\ = \SA\, hence \CB'\ = \CA\ 
due to symmetry and \CB'\ < \CB\ because B' > 120 c . Conversely, if \BC\ > \B'C\, 
then there exists B' G [B,S] with \CB'\ = \CA\, because \BC\ > \CA\ > \SC\. Then 
|AS| = \SB'\ < \SB\. 

Thus, the two-edges tree T = [A, B] U [B, C] is a minimal spanning tree for ABC, if 
and only if BC is the longest side of ABC, if and only if \AS\ < \BS\ and \AS\ < \CS\. 
Consider the points E G [B, S] and D G [C,S], such that \AS\ = \ES\ = \DS\, and 
put x = \AC\, y = \AB\, z = \DE\ = \AD\ = \AE\, and x' = \CD\, y' = \EB\. Then 
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Figure 18: To the calculation of sr3(R 2 , p 2 ). 



\SA\ = \SE\ = \SD\ = z/V3 and 

smt(M) = 3\SA\ + \DC\ + \EB\ = Viz + x' + y' and mst(M) = x + y, 

where M stands for the set {A, B, C}. But x < x' + z and y < y' + z, due to the triangle 
inequality, and hence 

i n r ^ V3z + x' +y' V3z + x' + y' V3z + x' + y'V3 

sr(M) = > = > . 

x + y x' + z + y' + z x' + y' + 2z 2 

Thus, we proved the following statement. 

Assertion 4.4 The following relation is valid: si3(M. 2 ,p 2 ) — v / 3/2. 

Remark. For small n it is already proved that sr„(R 2 , p 2 ) = (recently O. de Wet 

proved it for n < 7, see [16]). The proof of de Wet is based on the analysis of Du and 
Hwand method from |14j and understanding that it works for boundary sets with n < 7 
points. Also in 60th several lower bounds for sr(M 2 ,p2) were obtained, and the best of 
them is worse than in the third digit only. 

Problem 4.5 Very attractive problem is to prove that sr^R 2 ,^) = y/3/2, i.e. to prove 
Gilbert-Pollack Conjecture. The attempts to repair the proof of Du and Hwang have 
remained unsuccessful, so some fresh ideas are necessary here. 

4.3 Steiner Ratio of Other Euclidean Spaces and Riemannian 
Manifolds 

The following result is evident, but useful. 

Assertion 4.6 If Y is a subspace of a metric space X, i.e. the distance function on Y 
is the restriction of the distance function of X , then sr(Y) > sr(X). 

This implies, that sr(M",p 2 ) < sr(R 2 ,p 2 ) < Recall that Gilbert-Pollack con- 

jecture implies that the Steiner ratio of Euclidean plane attains at the vertex set of 
a regular triangle. In multidimensional case the situation is more complicated. The 
following result was obtained by Du and Smith [TT] 

Assertion 4.7 // M C 1" is the vertex set of a regular n-dimensional simplex, then 
sr(M) > sr(R", p 2 ) for n > 3. 
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Figure 19: Construction of the set P in R 3 (non-interesting but visual case n = 2). 



Proof. Consider the boundary set P in R n+1 , consisting of the following 1 + n(n + 1) 
points: one point (0, . . . , 0) and n(n + 1) points all whose coordinates except two are 
zero, one is equal to 1, and the remaining one is —1. It is clear that P is a subset of 
n-dimensional plane defined by the next linear condition: sum of all coordinates is equal 
to zero. Represent P as the union of the subsets P % — {x £ P | x 1 = 1} U {(0, . . . , 0)}. 
Notice that each set P l , i = 1, . . . ,n + 1, consists of n + 1 points and forms the vertex 
set of an regular n-dimensional simplex (to see that it suffices to verify that all the 
distances between the pairs of points from P l are the same and are equal to y/2). The 
configuration of 7 points in R 3 is shown in Figure [19] (this case is not important for us, 
but it is easy to draw). Now, mst(P) = (n + 1) mst(P z ), but for n > 3 we conclude that 
smt(P) < (n+1) smt(P z ), because the degree of the vertex (0, . . . , 0) in the corresponding 
network which is the union of the shortest networks for P % is equal to n + 1 > 4 that is 
impossible in the shortest network due to the Local Structure Theorem 13. II So, 

Taking as a heuristic for the length of a shortest tree connecting the vertex set of 
regular simplex the length of the network joining the center of the simplex with all its 
vertices we get the following estimate. 

Corollary 4.8 For any n > 3 the upper estimate 



is valid. 

One of the best general low estimates is obtained by Graham and Hwang in Q15] . 

Assertion 4.9 For any n>2 the lower estimate l/v3 < sr(R n ,p2) is valid. 

The best known upper estimate for R 3 is obtained by Smith and Smith [IB]. It is 
attained at an infinite boundary set which is known as "Smith sausage" and depicted 
in Figure [501 The corresponding value, obtained as the limit of the ratios for finite 
fragments, is as follows: 



283 3V2T 9^22- 2V2T 
700 ~ 700 + 140 ' 
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Figure 20: A finite fragment of infinite "Smith sausage". 



Notice that the idea of an infinite set is based on a deep result of Du and Smith estimating 
from below the number of points in a subset M of R™ such that sr(M) = sr(R™, p%) by a 
function f(n) rapidly increasing on n, see details in [17]. For example, /(50) = 53, but 
/(200) = 3 481 911. Therefore, it is difficult to expect to guess a finite set M in R™ with 
sr(M) = sr(R", p 2 ) for large n. 

Recently, the Steiner ratio of the Lobachevskii plane, and hence, of any Lobachevskii 
space has bin calculated by Innami and Kim, see [20] . 

Theorem 4.10 Steiner ratio of Lobachevskii space L n for any n > 2 is equal to 1/2. 

For general Riemannian manifold Ivanov, Cieslik and Tuzhilin, see |21| . obtained the 
following general result. 

Theorem 4.11 The Steiner ratio of n- dimensional Riemannian manifold is less than 
or equal to the Steiner ratio of the Euclidean space R" . 

5 Minimal Fillings 

This Section is devoted to minimal fillins, the third kind of optimal connections discussed 
in the Introduction. This problem appeared as a result of a synthesis of two classical 
problems: the Steiner problem on the shortest networks (it is discussed in Sections [3] 
and[4j, and Gromov's problem on minimal fillings. 

The concept of a minimal filling appeared in papers of Gromov, see [22]. Let M be a 
manifold endowed with a distance function p. Consider all possible films W spanning M, 
i.e., compact manifolds with the boundary M. Consider on W a distance function d that 
does not decrease the distances between points in M. Such a metric space W = (W, d) is 
called a filling of the metric space M. = (M, p) , see example in Figure [21] The Gromov 
Problem consists in calculating the infimum of the volumes of the fillings and describing 
the spaces W which this infimum is achieved at (such spaces are called minimal fillings). 

In the scope of Steiner problem, it is natural to consider M as a finite metric space. 
Then the possible fillings are metric spaces having the structure of one-dimensional strat- 
ified manifolds which can be considered as graphs whose edges have nonnegative weights. 
This leads to the following particular case of generalized Gromov problem. 

Let M be an arbitrary finite set, and G — (V. E) be a connected graph. We say, that 
G connects M or joins M, if M C V. Now, let A4 = (M, p) be a finite metric space, 
G = (V, E) be a connected graph joining M, and lj : E — > R + is a mapping into non- 
negative numbers, which is usually referred as a weight function and which generates the 
weighted graph Q = [G,ui). The function u generates on V the pseudo-metric d M (some 
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Figure 21: The space M is the circle S 1 with arc-metric. The films X in the both Figures 
are parts of the standard sphere containing M as a parallel. The left film X is not a 
filling since the distance between the points p and q in X is less than in M (the shortest 
path is depicted). The right film X is a filling of M. 

distances in a pseudo- metric can be equal to zero), namely, the ^-distance between the 
vertices of the graph Q is defined as the least possible weight of the paths in Q joining 
these vertices. If for any two points p and q from M the inequality p(p,q) < d LJ (p,q) 
holds, then the weighted graph Q is called a filling of the space M., and the graph G 
is referred as the type of this filing. The value mf(A4) = inf U)(Q), where the infimum 
is taken over all the fillings Q of the space M. is the weight of minimal filling, and each 
filling Q such that uj(Q) = mf(.M) is called a minimal filling. 

5.1 Parametric Networks and Optimal Connection Problems 

Here we give a common view on Steiner problem and minimal filling problem in terms 
of so-called parametric networks in a general metric space. 

Let X = (X, d) be a metric space and G = (V, E) be an arbitrary connected graph. 
Any mapping Y : V X is called a network in X parameterized by the graph G = (V, E), 
or a network of the type G. The vertices and edges of the network Y are the restrictions 
of the mapping Y onto the vertices and edges of the graph G, respectively. The length of 
the edge Y: vw — > X is the value d(T(v),T(w)j , and the length d(Y) of the network Y is 
the sum of lengths of all its edges. We shall consider various boundary value problems 
for graphs. To do that, we fix some subsets dG of the vertex sets V of our graphs 
G = (V, E), and we call such dG the boundaries. We always suppose that in each graph 
under consideration a boundary, possibly, an empty one, is chosen. The boundary dY of 
a network Y is the restriction of Y onto dG. If M C X is finite and M C r(V), then we 
say that the network Y joins or connects the set M. The vertices of graphs and networks 
which are not boundary ones are called interior vertices. The value 

smt(M) = inf {d(r) | Y is a network joining M} 

is called the length of shortest network for M. Notice that the network Y which joins 
M and satisfies d(Y) = smt(M) may not exist, see [53] and [23] for nontrivial examples. 
If such a network exists, it is called a shortest network connecting M, or for M . One 
variant of the Steiner problem is to describe the shortest networks for finite subsets of 
metric spaces H 

Now let us define minimal parametric networks in a metric space X — (X,d). Let 
G = (V,E) be a connected graph with some boundary dG, and let (p: dG — > X be a 

3 The denotation smt is an acronym for "Steiner Minimal Tree" which is a synonym for the shortest 
network whose edges are non-degenerate and, thus, it must be a tree. 
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mapping. By [G, ip] we denote the set of all networks F : V — > X of the type G such that 
dT = ip. We put 

mpn(G, tp) = inf d(T) 

re[G,<p] 

and we call this value the length of minimal parametric network. If there exists a network 
r 6 [G, if] such that d(T) = mpn(G, tp), then T is called a minimal parametric network 
of the type G with the boundary <p. 

Assertion 5.1 Let X = (X,d) be an arbitrary metric space and M be a finite subset of 
X . Then 

smt(M) = inf {mpn(G, cp) \ ip{dG) = M}, 

where the infimum is taken over all connected graphs G with a boundary dG and all 
mappings (p: dG —> X with ip(dG) = M . 

Thus, as in the case of the plane, the problem of calculating the length of the shortest 
network is reduced to investigation of minimal parametric networks. 

Let M. = (M, p) be a finite metric space and G = (V, E) be an arbitrary connected 
graph connecting M. In this case we always assume that the boundary of such G is fixed 
and equal to M. By Vl(M.,G) we denote the set of all weight functions u: E — > K such 
that (G, ui) is a filling of the space M.. We put 

mpf(M,G) = inf MG) 

and we call this value the weight of minimal parametric filling of the type G for the space 
A4. If there exists a weight function lu £ tt(A4,G) such that w(G) = mpf(A / i,G), then 
(G, uj) is called a minimal parametric filling of the type G for the space M.. 

Assertion 5.2 Let M = (M, p) be a finite metric space. Then 

mf(M) = inf{mpf(M,G)}, 

where the infimum is taken over all connected graphs G joining M . 

It is not difficult to show that to investigate shortest networks and minimal fillings 
one can restrict the consideration to trees such that all their vertices of degree 1 and 2 
belong to their boundaries. In what follows, we always assume that this condition 
holds, providing the opposite is not declared. 

To be more precise, we recall the following definition. We say that a tree is a binary 
one if the degrees of its vertices can be 1 or 3 only, and the boundary consists just of 
all the vertices of degree 1. Then each finite metric space has a binary minimal filling 
(possibly, with some degenerate edges), and a non-degenerate minimal filling (whose type 
is a tree and all whose vertices of degree 1 and 2 belong to its boundary in accordance 
with the above agreement), see [3]. 

5.2 Minimal Realization 

It turns out that the problem on minimal filling can be reduced to Steiner problem in 
special metric spaces and for special boundaries. 

Consider a finite set M = {p\, . . . ,_p n }, and let M = (M, p) be a metric space. We 
put pij = p(pi,pj). By we denote the n-dimensional arithmetic space with the norm 

||( W 1 ,..., U ")|| oo =max{| W 1 |,...,K|}, 

and by p^ the metric on generated by || • ||oo, i.e., poo{v,w) = \\w — v||oo. Let us 
define a mapping <pm ■ M —> RJ^ as follows: 

VMiPi) =Pi = (fin, ■ ■ ■ ,Pin)- 
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Assertion 5.3 The mapping ifM is an isometry with its image. 
Proof. This easily follows from the triangle inequality. Indeed, 

lift -ft || = max|p ifc - p jk \ > Pij, 

k 

because the value /9y stands at the ith and jth places of the vector ft — ft . On the other 
hand, pij > pik ~ Pjk for any k, due to the triangle inequality, hence ||ft — ft|| < pij, and 
Assertion is proved. 

The mapping ifM is called the Kuratowski isometry. 

Let Q = (G,(jj) be a filling of a space M. — (M, p), where G = (V, E), and be the 
pseudo-metric on V generated by the weight function uj. Denote by Em the edges set of 
the complete graph on M and put G — (V,E = E\J Em). Let Q be the weight function 
on E coinciding with metric p on Em and with u on E \ Em ■ Recall that do denotes the 
pseudo- metric on V generated by ui. 

We define the network Tg : V of the type G as follows: 

Tg(v) = (d (v,pi), . . .,du(v,p n )). 
This network is called the Kuratowski network for the filling Q . 
Assertion 5.4 We have dTg — (fiM- 

Proof. This easily follows from the filling definition. Indeed, the mapping dTg is defined 
on the set M only. By definition, 

r sOi) = {du,(Pi,Pl), ■ ■ ^doipi^n)), 

hence it suffices to show that d Q (pi,pk) = Pik for any k. The vertices pi and pu are joined 
by the edge PiPk of the weight pik in the graph G, and the weight of any other path in 
G connecting pi and pk is more than or equal to pik , because G is a filling. Assertion is 
proved. 

For any network T in a metric space (A, d) by ujt we denote the weight function on 
G induced by the network T, i.e., wr (vw) = d(T(v),T(w)). 

Corollary 5.5 Let Q = (G,cj) be a minimal parametric filling of a metric space (M,p) 
and T = Tg be the corresponding Kuratowski network. Then u) = LU-p . 

Let r be a network in a metric space X, let G be its parameterizing graph, and 
Ti. = (H, uj) be a weighted graph. We say that T and % are isometric, if there exists an 
isomorphism of the weighted graphs % and Q = (G,wr)- 

Corollary 15.51 and the existence of minimal parametric and shortest networks in a 
finite-dimensional normed space [T] imply the following result. 

Corollary 5.6 Let A4 = (Af, p) be a metric space consisting of n points, and lpm '■ M — ¥ 
be the Kuratowski isometry. For any graph G joining M there exists a minimal 
parametric filling of the type G of the space M.. Each minimal parametric filling of 
the type G of the space M. is isometric to the corresponding Kuratowski network, which 
is, in this case, a minimal parametric network of the type G with the boundary ^m- 
Conversely, each minimal parametric network of the type G on ifM (M) is isometric to 
some minimal parametric filling of the type G of the space M. . 

Corollary 5.7 Let Ai = (M, p) be a metric space consisting of n points, and tp^ : M — > 
R.^ be the Kuratowski isometry. Then there exists a minimal filling Q for M., and the 
corresponding Kuratowski network Tg is a shortest network in the space joining the 
set ipm(M). Conversely, each shortest network on Lpj^{M) is isometric to some minimal 
filling of the space Ai . 
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5.3 Minimal Parametric Fillings and Linear Programming 

Let M = (M,p) be a finite metric space connected by a (connected) graph G = (V,E). 
As above, by Q(A4, G) we denote the set consisting of all the weight functions oj : E — >• R + 
such that Q = (G. lu) is a filling of M., and by £l m (A4, G) we denote its subset consisting 
of the weight functions such that Q is a minimal parametric filling of M. . 

Assertion 5.8 The set ft(A4,G) is closed and convex in the linear space M. E of all the 
functions on E, and tt m (A4,G) C tt(M,G) is a nonempty convex compact. 

Proof. It is easy to see, that the set ft(A4,G) C M. E is determined by the linear in- 
equalities of two types: w(e) > 0, e <E E, and ~^2 ef z lpq w(e) > p(p,q), where j pq stands 
for the unique path in the tree G connecting the boundary vertices p and q. Therefore, 
£l(A4,G) is a convex closed polyhedral subset of R that is equal to the intersection 
of the corresponding closed half-spaces. The weight functions of minimal parametric 
fillings correspond to minima points of the linear function X) e e-E w ( e ) restricted to the 
set Q,{Ai,G). Thus, the problem of minimal parametric filling finding is a linear pro- 
gramming problem, and the set £l m (A4 , G) of all minima points is a nonempty convex 
compact polyhedron (the boundedness and, hence, compactness of this set follows from 
increasing of the objective function with respect to each its variable). 

5.4 Generalized Fillings 

Investigating the fillings of metric spaces, it turns out to be convenient to expand the 
class of weighted trees under consideration permitting arbitrary weights of the edges (not 
only non- negative) . The corresponding objects are called generalized fillings, minimal 
generalized fillings and minimal parametric generalized fillings. Their weights for a metric 
space M and a tree G are denoted by mf _ (M ) and mpf _ (M , G) , respectively. 

For any finite metric space A4 — (M, p) and a tree G connecting M, the next evident 
inequality is valid: mpf_(A^,G) < xapi(M.,G). And it is not difficult to construct 
an example, when this inequality becomes strict, see Figure [2"2l However, for minimal 
generalized fillings the following result holds, see [2"5] . 

Theorem 5.9 (Ivanov, Ovsyannikov, Strelkova, Tuzhilin) For an arbitrary finite 
metric space M. , the set of all its minimal generalized fillings contains its minimal filling, 
i.e. a generalized minimal filling with nonnegative weight function. Hence, mf_(A4) = 
mf(M). 

5.5 Formula for the Weight of Minimal Filling 

Let M. = (M, p) be a finite metric space, and G be a tree connecting M. Choose an 
arbitrary embedding G' of the tree G into the plane. Consider a walk around the tree 
G' . We draw the points of M consecutive with respect to this walk as a consecutive 
points of the circle S 1 . Notice that each vertex p from M appears degp times. For 
each vertex p G M of degree more than 1, we choose just one arbitrary point from the 
corresponding points of the circle. So, we construct an injection v: M — > S . Define 
a cyclic permutation 7r as follows: 7r(p) = q, where v{q) follows after v{p) on the circle 
S . We say that ir is generated by the embedding G' (this procedure is not unique due to 
different possible choices of v). Each 7r generated in this manner is called a tour of M 
with respect to G. The set of all tours on M with respect to G is denoted by 0(M, G). 
For each tour n e 0(M, G) we put 

p(M,G,ir) = 2 p ( x > 7r ( x )) 
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Figure 22: Minimal parametric filling (left) and minimal generalized parametric filling 
(right) of the vertex set of the plane rectangle with sides 3 and 4. The type is the same: 
the moustaches connects the diagonal pairs of the vertices. The interior edge has to be 
zero in the case of the filling and can be negative in the case of the generalized filling. 
Here 9 = mpf_(M, G) < mpf(.M, G) = fO. 

and we call this value by the half-perimeter of the space M. with respect to the tour 7r. 
The minimal value of p(Ai,G, tt) over all tt £ 0(M,G) for all possible G (in fact, over 
all possible cyclic permutations tt on M) is called the half-perimeter of the space A4. 
A. Ivanov and A. Tuzhilin proposed the following hypothesis. 

Conjecture 5.10 For an arbitrary metric space Ai — {M,p) the following formula is 
valid 

mi(A4) = min max p(AA,G, tt), 

G -rr£0{M,G) 

where minimum is taken over all binary trees G connecting M . 

A. Yu. Eremin 26 constructed a counter-example to the Conjecture 15 . 101 and showed 
that if one changes the concept of tour by the one of multitour, introduced by him, then 
the Conjecture 15. 101 holds. 

To define the multitours, let us consider the graph in which every edge of G is taken 
with the multiplicity 2k, k > 1. The resulting graph possesses an Euler cycle consisting 
of irreducible boundary paths — the ones which do not contain properly other boundary 
paths. This Euler cycle generates a bijection tt: X — >• X, where X = U^ =1 M, which is 
called multitour of M with respect to G, see an example in Figure 1231 The set of all 
multitours on M with respect to G is denoted by C M (M, G). 

Let M. = (M, p) be a finite metric space, and G be a tree connecting M. As in the 
case of tours, for each multitour tt £ O^iM, G) we put 

p(M,G, k) = ^J2 p( x > n W) ■ 

Theorem 5.11 (A. Yu. Eremin) For an arbitrary finite metric space Ai — (M,p) and 
an arbitrary tree G joining M , the weight of minimal parametric generalized filling can 
be calculated as follows 

mpf_ (M,G) = max{p(X,G,vr) | tt £ O ^ [M , G) } . 

The weight of minimal filling can be calculated as follows 

m£(M) = mf_(M) = minmax{p(7W,G,7r) | tt £ 0^(M,G)}, 

G 

where minimum is taken over all binary trees G connecting M . 
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Figure 23: A part of a moultitour with multiplicity 2 (left), and the irreducible boundary 
paths the multitour consists from (right). The multitour starts as a green polygonal line 
and becomes blue when multiplicity of edges becomes more than 2. 

5.6 Minimal Fillings for Generic Metric Spaces 

Theorem 15.111 gives an opportunity to get several interesting corollaries. To formulate 
one of them, we need to define what is a "generic" metric space. Notice that the set of 
all metric spaces consisting of n points can be naturally identified with a convex cone in 
fljn(n-i)/2 suffices to enumerate the set of all two-elements subsets of these spaces and 
assign to each such space the vector of the distances between the pairs of points). This 
representation gives us an opportunity to speak about topological properties of families 
of metric spaces consisting of a fixed number of points. 

We say, that some property holds for a generic metric space, if for any n this property 
is valid for an everywhere dense set of n-point metric spaces. 

The following result can be found in 2(> . 

Corollary 5.12 (A. Yu. Eremin) Each general finite metric space has a minimal fill- 
ing which is a nondegenerate binary tree. 

5.7 Additive Spaces and Minimal Fillings 

The additive spaces are very popular in bioinformatics, playing an important role in 
evolution theory and, more general, in an hierarchy modeling. Recall that a finite metric 
space A4 = (M, p) is called additive, if M can be joined by a weighted tree Q — (G, w) 
such that p coincides with the restriction of du onto M. The tree Q in this case is called 
a generating tree for the space M.. 

Not any metric space is additive. An additivity criterion can be stated in terms of 
so-called 4 points rule: for any four points pi, pj, pk, pi, the values p(pi,Pj) + p(pk,Pi), 
p{Pi,Pk) + p(PjtPl)t piPiiPi) + p(Pj,Pk) are the lengths of sides of an isosceles triangle 
whose base does not exceed its other sides. 

Theorem 5.13 (|28j, [29j , |30j, |31j) A metric space is additive, if and only if it sat- 
isfies the 4 points rule. In the class of non-degenerate weighted trees, the generating tree 
of an additive metric space is unique. 

The next criterion solves completely the minimal filling problem for additive metric 
spaces. 

Theorem 5.14 Minimal fillings of an additive metric space are exactly its generating 
trees. 
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The next additivity criterion is obtained by O. V. Rubleva, a student of Mechanical 
and Mathematical Faculty of Moscow State University, see [32] . 

Assertion 5.15 (O. V. Rubleva) The weight of a minimal filling of a finite metric 
space is equal to the half-perimeter of this space, if and only if this space is additive. 

In the scope of Assertion l5.15[ we conjectured that if there exists a tree connecting a 
metric space such that all the corresponding half-perimeters are equal to each other, then 
the space is additive. It turns out that it is not true. Z. N. Ovsyannikov suggested to 
consider a wider class of spaces, so called pseudo-additive spaces, for which our conjecture 
becomes true, see [27] . 

A finite metric space Ai = (M, p) is said to be pseudo-additive, if the metric p 
coincides with d^ for a generalized weighted tree (G,uj) (which is also called generat- 
ing), where the weight function oj can take arbitrary (not necessary nonnegative) val- 
ues. Z.N. Ovsyannikov shows that these spaces can be described in terms of so-called 
weak 4-points rule: for any four points pi, pj, pk, pi, the values p(pi,Pj) + p(pk,Pi), 
p{pi,pu) + p{Pj,Pi), p{Pi,Pi) + p{Pj,Pk) are the lengths of sides of an isosceles triangle. 
The generating tree is also unique in the class of non-degenerate trees. Moreover, the 
following result is valid, see [27j . 

Theorem 5.16 (Z.N. Ovsyannikov) Let M. = {M,p) be a finite metric space. Then 
the following statements are equivalent. 

• There exist a tree G such that M coincides with the set of degree 1 vertices of G 
and all the half-perimeters p{M, G, 7r) of M corresponding to the tours around G 
are equal to each other. 

• The space M. is pseudo-additive. 

Moreover, the three G in this case is a generating tree for the space M. . 

It would be interesting to see, what role these pseudo-additive spaces could play in 
applications. 



5.8 Examples of Minimal Fillings 

Now let us give several examples of minimal filling and demonstrate how to use the 
technique elaborated above. 



5.8.1 Triangle 

Let M. = (M,p) consist of three points p\, P2, and p%. Put pij = p(pi,pj). Consider the 
tree G = (V,E) with V = M U {v} and E = Define the weight function uj on 

E by the following formula: 

, ■. _ Pij + Pik - Pjk 

where {i,j, k} — {1, 2, 3}. Notice that d^ restricted onto M coincides with p. Therefore, 
M. is an additive space, Q = (G, w) is a generating tree for A4, and, due to Theorem l5.14[ 
Q is a minimal filling of A4 . 

Recall that the value (pij + pik — Pjk)/2 is called by the Gromov product (pj,Pk)pi of 
the points pj and pk of the space A4 with respect to the point pi, see |33j . 
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5.8.2 Regular Simplex 



Let all the distances in the metric space M. are the same and are equal to d, i.e. M. is 
a regular simplex. Then the weighted tree Q — (G,uj), G = (V,E), with the vertex set 
V = M U {v} and edges vm, m £ M, of the weight d/2 is generating for M. Therefore, 
the space M. is additive, and, due to Theorem l5.14l Q is its unique nondegenerate minimal 
filling. If n is the number of points in M, then the weight of the minimal filling is equal 
to dn/2. 

5.8.3 Star 

If a minimal filling Q = (G, u) of a space M. = (M, p) is a star whose single interior 
vertex v is joined with each point £ M, 1 < i < n, n > 3, then the metric space M. 
is additive [3]. In this case the weights of edges can be calculated easily. Indeed, put 
e s ; = vpi. Since a subspace of an additive space is additive itself, then we can use the 
results for three-points additive space, see above. So, we have w(ej) = (pj,Pk)pi-> where 
Pi, pj, and pk are arbitrary distinct boundary vertices. 

5.8.4 Mustaches of Degree more than 2 

Let G — (V, E) be an arbitrary tree, and v £ V be an interior vertex of degree (fe + 1) > 3 
adjacent with k vertices u>i, . . . , Wk from dG. Then the set of the vertices {w\, . . . ,Wk}, 
and also the set of the edges {vwi, . . . , vwk}, are referred as mustaches. The number k is 
called by the degree, and the vertex v is called by the common vertex of the mustaches. 
An edge incident to v and not belonging to {vw\, . . . , vwk} is called the root edge of the 
mustaches under consideration. 

As it is shown in [3] , any mustaches of a minimal filling of a metric space forms an 
additive subspace. If the degree of such mustaches is more than 2, then we can calculate 
the weights of all the edges containing in the mustaches just in the same way as in the 
case of a star. 

5.8.5 Four-Points Spaces 

Here we give a complete description of minimal fillings for four-points spaces, see details 

»> m 

Proposition 5.17 Let M = {pi,p2,P3,P4}, and p be an arbitrary metric on M. Put 
Pij = p(pi,pj). Then the weight of a minimal filling Q — (G,w) of the space M. = {M,p) 
is given by the following formula 

i(min{p 12 + p 34 , P13 + P24, Pu + P23} + max{p 12 + p 3 4, P13 + P24, Pu + £23})- 

If the minimum in this formula is equal to p^ + p rs , then the type of minimal filling is 
the binary tree with the mustaches {pi,pj} and {p r ,p s }. 

We apply the obtained result to the vertex set of a planar convex quadrangle. 

Corollary 5.18 Let M be the vertex set of a convex quadrangle P1P2P3P4 C K 2 and 
p(Pi>Pj) = \\Pi~Pj\\- The weight of a minimal filling of the space {M,p) is equal to 



The topology of minimal filling is a binary tree with mustaches corresponding to opposite 
sides of the less total length. 




29 



5.9 Ratios 



The Steiner ratio is discussed in Section 01 Here we define some other ratios based on 
minimal fillings, which could be more available for calculating, and which could be useful 
to calculate the Steiner ratio, as we hope. 

5.10 Steiner— Gromov Ratio 

For convenience, the sets consisting of more than a single point are referred as nontrivial. 
Let X = (X,p) be an arbitrary metric space, and let M C X be some finite subset. 
Recall that by mst(M, p) we denote the length of minimal spanning tree of the space 
(M, p). Further, for nontrivial M, we define the value 



and call it the Steiner- Gromov ratio of the subset M. The value infsgr(Af), where the 
infimum is taken over all nontrivial finite subsets of X, consisting of at most n vertices is 
denoted by sgv n (X) and is called the degree n Steiner- Gromov ratio of the space X. At 
last, the value inf sgr n (<f), where the infimum is taken over all positive integers n > 1 is 
called the Steiner- Gromov ratio of the space X and is denoted by sgr(X), or by sgr(X), 
if it is clear what particular metric on X is considered. Notice that sgr n (X) is a non- 
increasing function on n. Besides, it is easy to see that sgv 2 {X) = 1 and sgr 3 (A?) =3/4 
for any nontrivial metric space X. 

Assertion 5.19 The Steiner-Gromov ratio of an arbitrary metric space is not less than 
1/2. There exist metric spaces whose Steiner-Gromov ratio equals to 1/2. 

Recently, A. Pakhomova, a student of Mechanical and Mathematical Faculty of 
Moscow State University, obtained an exact general estimate for the degree n Steiner- 
Gromov ratio, see |35j . 




Assertion 5.20 (A. Pakhomova) For any metric space X the estimate 



is valid. Moreover, this estimate is exact, i.e. for any n > 3 there exists a metric space 
X n such that sgv n (X n ) = n/(2n — 2). 

Also recently, Z. Ovsyannikov [36] investigated the metric space of all compact subsets 
of Euclidean plane endowed with Hausdorff metric. 

Assertion 5.21 (Z. Ovsyannikov) The Steiner ratio and the Steiner-Gromov ratio of 
the metric space of all compact subsets of Euclidean plane endowed with Hausdorff metric 
are equal to 1/2. 

5.11 Steiner Subratio 

Let X = (X,p) be an arbitrary metric space, and let M C X be some its finite subset. 
Recall that by smt(M, p) we denote the length of Steiner minimal tree joining M. Further, 
for nontrivial subsets M, we define the value 



and call it by the Steiner subratio of the set M. The value inf ssr(M), where infimum is 
taken over all nontrivial finite subsets of X consisting of at most n > 1 points, is denoted 



sgr(M) 



mf(M, p)/mst(M, p) 




ssr(M) 



mf(M, p)/smt(M, p) 
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by ssr„(A') and is called the degree n Steiner subratio of the space X. At last, the value 
infssr„(A'), where the infimum is taken over all positive integers n > 1, is called the 
Steiner subratio of the space X and is denoted by ssr(A"), or by ssr(X), if it is clear what 
particular metric on X is considered. Notice that ssr n (X) is a nonincreasing function on 
n. Besides, it is easy to see that ss^A") = 1 for any nontrivial metric space X. 

Proposition 5.22 ssr 3 (tt") = y/3/2. 

The next result is obtained by E. I. Filonenko, a student of Mechanical and Mathe- 
matical Department of Moscow State University, see [34] . 

Proposition 5.23 (E.I. Filonenko) ssr 4 (IR 2 ) = V3/2. 

Conjecture 5.24 The Steiner subratio of the Euclidean plane is achieved at the regular 
triangle and, hence, is equal to \/2>j1. 

Recently, A. Pakhomova obtained an exact general estimate foe the degree n Steiner 
subratio, see [35] . 

Proposition 5.25 (A. Pakhomova) For any metric space {X} the estimate 

ssr n (X) > 

v 1 ~ In- 2 

is valid. Moreover, this estimate is exact, i.e. for any n > 3 there exists a metric space 
X n such that ssr n (A? n ) = n/(2n — 2). 

Also recently, Z. Ovsyannikov [36] investigated the metric space of all compact subsets 
of Euclidean plane endowed with Hausdorff metric. 

Proposition 5.26 (Z. Ovsyannikov) Let C be the metric space of all compact subsets 
of Euclidean plane endowed with Hausdorff metric. Then ssr3(C) = 3/4 and ssr4(C) = 
2/3. 
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